Interference between independent fluctuating condensates 
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We consider a problem of interference between two independent condensates, which lack true 
long range order. We show that their interference pattern contains information about correlation 
functions within each condensate. As an example we analyze the interference between a pair of 
one dimensional interacting Bose liquids. We find universal scaling of the average fringe contrast 
with system size and temperature that depends only on the Luttinger parameter. Moreover the 
full distribution of the fringe contrast, which is also equivalent to the full counting statistics of the 
interfering atoms, changes with interaction strength and lends information on high order correla- 
tion functions. We also demonstrate that the interference between two-dimensional condensates at 
finite temperature can be used as a direct probe of the Kosterlitz-Thouless transition. Finally, we 
discuss generalization of our results to describe the intereference of a periodic array of independent 
fluctuating condensates. 



I. INTRODUCTION 

An important property of Bose Einstein condensates 
(BEC) is the existence of a coherent macroscopic phase. 
Thus, a crucial benchmark in the study of such sys- 
tems was the observation of interference fringes when 
two independent condensates were let to expand and 
overlap^. This "two slit" experiment was carried out with 
cold atoms in three dimensional harmonic traps, where a 
true condensate exists. The interference fringe amplitude 
should then be proportional to the condensate fraction, 
as was indeed observed. However, with current trapping 
technology it is possible to confine the bosonic atoms 
to oneSiii or two dimensions^ where a true condensate 
may not exist. Instead, these systems are characterized 
by off-diagonal correlations that decay as a power-law 
or exponentially in space. What is the interference pat- 
tern that arises when two such imperfect condensates are 
let to expand and overlap? This question is not just of 
general academic interest. Recently there have been a 
number of beautiful experiments where the interference 
between independent co nde nsates was directly observed, 
see for example Refs. 0[7ll3l- 

Here we address this problem theoretically and show 
that the result depends crucially on the correlations 
within each condensate. Therefore, such an experiment 
would provide a direct and simple probe of the spatial 
phase correlations. In principle spatial phase correlations 
may also be extracted from juggling experimentsSiifliii or 
the momentum distribution measured by the free expan- 
sion of a single condensate^. However, creating strongly 
interacting low dimensional systems typically requires us- 
ing low density atomic gases, which makes juggling ex- 
periments very challenging. Also the highly anisotropic 
expansion of low dimensional condensates inhibit mea- 
surements of the momentum distribution in the slowly 
expanding longitudinal direction. A method for probing 
the phase correlations directly in real space would there- 
fore be very useful. 



II. RESULTS AND DISSCUSSION 

The simplest geometry that we consider is illustrated 
in Fig. ^ It consists of two parallel one dimensional 
(ID) condensates a distance d apart. After the atoms 
are released from the trap they are let to expand to a 
transverse size much larger than d, while no significant 
expansion occurs in the axial direction. An absorption 




FIG. 1: Schematic view of the possible experimental setup, 
which produces interference pattern between two independent 
ID condensates. Here L and d are the imaging length and the 
separation between the condensates, respectively. 

image is then taken by a probe beam directed along the 
condensate axis. A similar setup is considered for two di- 
mensional (2D) condensates on parallel planes (see Fig. 
|2Jl. As usual, the absorption image gives the instanta- 
neous three dimensional density profile, integrated along 
the beam axis. p{x) = dz al^j{x, z) atof{x, z), where 

al^j are the Bose creation operators with the subscript 
"to/" , emphasizing that the corresponding operators are 
taken after free expansion of atoms, z is the axial co- 
ordinate and X is the coordinate along the detector (see 
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ing component at wavector Q = md/Ht 



{p{xi)p{x2))int = (|Aq|^ 




(2) 



FIG. 2: Same as in FigQbut for two dimensional condensates. 



Fig.nj. The length L is typically given by the focal length 
of the imaging beam. It may also be controlled more pre- 
cisely by applying magnetic field gradients, so that only 
a specified section of the cloud is resonant with the probe 
light. In principle, one can consider an experiment with 
a probe beam orthogonal to the plane containing two 
parallel ID condensates. In this case, it is possible to 
integrate the resulting interference image within an arbi- 
trary interval and obtain dependence of the interference 
contrast on L (note that this dependence characterizes a 
single run and a series of experiments is still needed to 
find the average contrast). The other advantage of this 
set-up is that it can reveal the presence of dipolar oscilla- 
tions in individual condensates. These modes correspond 
to the center of mass motion and are not affected by in- 
teractions. Dipolar oscillations induce an overall tilt in 
the interference peak position and can be easily removed 
by integrating p{x) along a line tilted with respect to the 
z axis. However, since most of the current experimental 
systems do not allow imaging beams that are perpen- 
dicular to ID condensates, we concentrate on the setup 
shown in Fig. ^ 

To discuss the interference contrast we consider the 
correlation function of the density operator 



{p{xi)p{x2)) = S{xi - X2) / dzp{z) 



dzidz2 {al^f in )al^j: ir2)atof {n )atof (^2 )) , ( 1) 



where stands for (xi,Zi). Single particle oper- 
ators in should be taken after the expansion 
time t. We can relate them to operators before 
the expansionji^ atof{x,z) — ai{z)e^^^''^^^~^'^^*^^™ + 
a2(2;)e*'^^'^^)^~*'^i*/^'", with ai_2 being operators in the 
two condensates and Q12 = m{x ± d/2)/ht. We thus 
find that the correlation function in Q has an oscillat- 



dzidz2{al{zi)ai{zi)a\{z2)a2{z2))- (3) 



Here Aq = J dza\{z)a2{z) is the quantum observ- 
able corresponding to the amplitude of the interference 
fringes. It can be extracted from the time of flight (TOF) 
absorption image by taking the Fourier transform of the 
density profile. Alternatively one can directly probe Aq 
by studying the oscillating component in the density au- 
tocorrelation function. Both methods were successfully 
used in recent experimentsi^ii^. In practice it might be 
easier to study the interference contrast rather than the 
absolute value of the fringe amplitude. In this case one 
has to divide Aq by the imaging length L. If the two con- 
densates are decoupled from each other, the expectation 
value of (Aq) vanishes. This docs not mean that \Aq\ 
is zero in each individual measurement but rather shows 
that the phase of Aq is randomi^. Said differently, Aq 
is finite in each experimental run, but its average over 
many experiments vanishes. To determine the ampli- 
tude of interference fringes in individual measurements 
one should consider an expectation value of the quantity 
that does not involve the random phase of Aq . This nat- 
urally brings us to equation ©. From shot to shot |AqP 
fluctuates as well, and Eq. gives its average value. 

If the two condensates are identical (but still indepen- 
dent) we may simplify expression 



{\Aq\ 



L 



dz{aHz)a{0)y 



(4) 



Here we neglected boundary effects by integrating over 
the center of mass coordinate and assuming that the cor- 
relations depend only on (zi — Z2). Eq. |0J can be gener- 
alized for the case of parallel 2D condensates by taking 
z to represent the planar coordinates. 

To gain intuition into the physical meaning of the 
fringe amplitude let us first address two simple limiting 
cases. First, consider the situation where {a^{z)a{0)) de- 
cays exponentially with distance with a correlation length 
^ << L. Then Eq. (gj) implies that \Aq\ cx ^/L^, 
which has a very simple physical meaning. Since the 
phase is only coherent over a length ^ the system is ef- 
fectively equivalent to parallel chains with L/^ pairs of 
independent condensates. Each pair contributes inter- 
ference fringes with a constant amplitude proportional 
to ^ and a random phase. The total amplitude Aq is 
therefore the result of adding L/S^ independent vectors 
of constant length ^ and random direction, hence we get 
^/L^ scaling. Note that the interference contrast, which 
is proportional to Ag/L, the ratio of fringe amplitude to 
the background signal, scales as -^/^/L. This observation 
is similar in spirit to that made in Ref. Q of interference 
between 30 independent condensates in a chain. Fringes 
can be seen, though their average amplitude is suppressed 
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by a factor of vSO compared to the interference between 
two condensates. Now consider the opposite hmit of per- 
fect condensates, for which (at(z)a(O)) is constant. In 
this case Eq. Q imphes that \Aq\ (x L. PictoriaUy this 
is the result of adding vectors with a uniform phase, re- 
suhing in a fringe amphtude which scales as the total 
size of the system. This is essentially the result of the 
experiment in Ref. jlj. 



A. One dimensional Bose liquids 

We proceed to discuss the case of a ID interacting gas. 
We first consider a system at sufficiently low temper- 
ature, >> L, where is the temperature depen- 
dent correlation length defined in Eqs. H17fl and H18() be- 
low. In this regime the correlations decay as a power law 
rather than exponentially. We therefore, expect that the 
fringe amplitude will somehow interpolate between the 
two simple limits considered above. Specifically, at long 
wavelengths the ID Bose gas is described by a Luttinger 
liquid^^: and the long distance off-diagonal correlations 
behave as 

(at(z)a(0))~p(^^j . (5) 

Here p is the particle density, is the healing length, 
which also serves as the short range cutoff, and K is the 
Luttinger parameter. For bosons with a repulsive short- 
range potential, K ranges between 1 and oo, with K = 1 
corresponding to strong interactions, or "impenetrable" 
bosons while if — > oo for non interacting bosons. Sub- 
stituting Eq. (O into Eq. Q and assuming that L ^ S^fi 
we arrive at one of our main results 



(6) 



where C is a constant of order unity. Thus we see that the 
amplitude of the interference fringes {Aq = a/ {\Aq p)) , 
scales with a non trivial power of the imaging length. In 
the non interacting limit {K — > oo) the scaling is linear 
Aq ^ L as expected for a fully coherent system. In- 
terestingly Aq ~ -s/i in the hard core limit {K — 1), 
as in systems with short range correlations which were 
discussed above. A more careful examination of the inte- 
gral shows that at if = 1 there are additional logarithmic 
corrections to the power law scaling 

Having derived the amplitude of interference fringes, 
an interesting question is how this amplitude fluctu- 
ates from one experimental run to the next. To an- 
swer this question one should consider higher moments 
of the operator l^gp. We find that all moments have 
the general form: (|Aq|2") = {\Aq\^)^T,,{K). It fol- 
lows that when |AqP is normalized, its distribution func- 
tion P(|Aq|V(|^qP)) is fully determined by the Lut- 
tinger parameter K. In particular for large K the func- 
tion P(|Aq|V(|^qH) becomes very narrow character- 

-(|AqP)2/(|AqP) « 



7r/(\/6if). As expected, for if — > cx) the width of the 
distribution goes to zeroii. In the opposite limit of 
if — !■ 1 the distribution takes another simple limiting 
form: P{x) « e~^. It is interesting to point out that 
not only the scaling with L but the whole distribution 
of the interference fringes at if 1 is equivalent to the 
one that arises in systems with short range correlations. 
In Fig. Owe plot schematically the two distributions ob- 
tained in the opposite limits of strong and weak inter- 
actions. Note that the true distribution at large K is 
slightly asymmetric. We address this issue in a separate 
paper^Sj. We see that there is an interesting crossover 
from a narrow distribution in weakly interacting systems 
[K 3> 1) to a wide poissonian distribution in the hard- 
core limit (if — > 1). In fact, one can show that there is a 
close connection between the distribution function P(x) 
and the partition function of the boundary Sine-Gordon 
problem. This gives a direct link between interference 
of two ID quasi-condensates and thermodynamics of a 
single impurity problem in a Luttinger liquid. We defer 
the details of this problem and the crossover from weak 
to strong interactions for future investigation>i^. 




ized by the width a = ^ ^Aq 



FIG. 3: Schematic plots of the distribution function P(x) 
of the normalized intensity of interference fringes [x — 
|AqP/(|AqP)) in strongly and weakly interacting regimes: 
if ^ 1 and if = 10. 



We showed above that it is possible to extract if by 
analyzing either the scaling of the fringe amplitude with 
system size or by analyzing its distribution at a given sys- 
tem size. Another approach involves changing the angle 
9 between the probe beam and the condensate axis while 
keeping the imaging length fixed. The resulting absorp- 
tion image then corresponds to integration of the cloud 
density along a line at an angle to the z axis. Then 
Eq. should be changed accordingly. Let x be the 
transverse coordinate on the screen (z = 0), then in the 
second term in the right hand side of ^ we substitute 
X = X — ztanO. On squaring we obtain the analogue of 
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Eq. Q 

{\Aq\^{9)) =L ( dz{a\z)am^cos[q{e)z], (7) 

where q{6) = fcotan^. For sufficiently large imaging 
length {qL >> 1) this yields 

{\AQ\\e))^{p^Lin){iHqf"'-\ (8) 

Thus the Luttinger parameter may be extracted from the 
angle dependence of the fringe amplitude. For qL << 
1 (i.e. very small angles) Eq. (O reduces to Eq. 
Note that if one uses the imaging beam orthogonal to the 
condensates, then 9 will be simply the angle between the 
z axis and the direction of integration of the interference 
contrast. 

Before concluding this section let us address the ef- 
fect of temperature. It is well known that at any finite 
temperature the off-diagonal correlations in a one dimen- 
sional Bose system must be short ranged. Specifically at 
sufficiently long distances they decay exponentially with 
a correlation length ~ l/T. The zero temperature 
results presented above are valid at sufficiently low tem- 
perature that >> L. At higher temperature such that 
<< L, the scaling of the fringe amplitude with length 
must be \Aq\ ~ \/L as discussed above for systems with 
short range correlations. However, as long as >> S,h 
the temperature dependence of the fringe amplitude is 
universal and it depends only on the Luttinger parame- 
ter: 

(l^«r-V6.(|.^)"*- W 

This provides another experimental method to extract 
Luttinger physics. 

B. Two dimensional systems 

We now consider a pair of parallel 2D condensates. In 
direct analogy to the ID condensates, the imaging axis 
may be taken parallel or at some angle to the plane of 
the condensates. In the former case one should consider 
the scaling of the fringe amplitude with imaging length, 
while in the latter case the variation with angle. 

It is well known that in two dimensions, long range 
order may exist only at zero temperature. At sufficiently 
low temperatures off-diagonal correlations are algebraic, 
with 

{aHr)am ^ p(^^y . (10) 

for r >> ^/j. On the other hand, above the Kosterlitz- 
Thouless (KT) transition aX T = Tc, the correlations 
decay exponentially. We will show that this transition 
is characterized by a jump in the behavior of the fringe 



amplitude, related to the well known universal jump of 
the superfluid stiffness at Tc- 

The exponent in (|l()|l is given by a = mT /2-Kps{T)h^ . 
For weakly interacting bosons at temperatures well be- 
low Tc, Ps{T) is simply equal to the density p. As one 
approaches the transition pa is renormalized by fluctua- 
tions, and at the transition Ps{Tc) = 'imTc/'KfT? . There- 
fore the exponent a assumes a universal value ac = 1/4 
at the transition. Thus for temperatures T < Tc we have 
< a < ac- 

Let us now discuss the consequences of this physics to 
the experimentally measurable fringe amplitude. As il- 
lustrated in Fig. 12 the interference pattern is now truly 
2D in the sense that cuts along x at different coordinate 
y display a different fringe pattern. To obtain a ID pat- 
tern as a function of x alone we may integrate the image 
intensity over an "integration length" Ly. Recall that 
in addition the imaging process automatically integrates 
over an imaging length along the z axis. Now the gen- 
eralization of Eq. Q to the 2D case is straight forward 

Ly 

{\Aq\^) ^ LyL, J dy J dz{a^{y,z),a{0,0))\ (11) 



For simplicity we assume that Ly and L^ are scaled si- 
multaneously as Ly ~ Lz — V^, with fl the imaging 
area. Then using Eqs. pO(l and 1111) we find that for 
T<Tc 

{\AQf)^n'-". (12) 

So below the transition the scaling of \Aq\ with size 
ranges from linear at T = to fi^-*^^ at the KT transition 
(ckc = 1/4). On the other hand, for T > Tc the corre- 
lations decay exponentially and \Aq\ ^ \/f2. Hence we 
find a universal jump at T = in the power character- 
izing the size dependence of \Aq\. One can also consider 
a setup where only dependence on one length, say Ly, is 
studied while the other one Lz is fixed. Then if Ly^L^, 
Eq. ifT^ immediately generahzes to (|v4qP) ~ ^y^^"- In 
this case the power jump in \Aq\ is bigger: from 0.75 to 
0.5 as T crosses Tc- This jump is a direct signature of 
the KT physics, to be contrasted with the result for ID 
condensates, where the scaling power with system size 
interpolated smoothly between 1 and 2. It should be 
noted, however, that a ID system on an optical lattice, 
which undergoes a Mott transition at T = 0, would dis- 
play a universal jump similar to the 2D case discussed 
here. In the same way one can study the shape of the 
distribution of the interference amplitude and find that 
as T increases to Tc the distribution gradually broadens 
but always remains relatively narrow. On the other hand, 
as T becomes larger than Tc the distribution assumes a 
broad poissonian form. 

The analysis for imaging the 2D condensates with a 
slanted probe beam can be carried over from the ID case. 
The scaling of the interference contrast with q = ko tan 9, 
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at constant imaging area, is then (|^qP) ~ 1/q^ ^" be- 
low the KT transition, and {\Aq\'^) ~ 1/(1 + q^i'^f'^ 
above it. Again the transition is characterized by a uni- 
versal jump of the power at small q. We emphasize that 
9 can be either the angle between the beam and the z 
axis (see Fig. [JJ or the angle between the y axis and the 
direction of integration. The latter is preferable, because 
within a single experimental shot it is possible to obtain 
the whole angular dependence of Ag. 

Regardless of the experimental approach of choice, the 
interference between a parallel pair of independent 2D 
condensates can serve as a direct probe of Kosterlitz- 
Thouless physics. However a word of caution is in order. 
The correlation length, which coincides with the healing 
length at very low temperaturesiS (^t ~ S,h = ^1 \J'nigp), 
diverges at the KT transition as cx exp(&/-v/rc — T)^. 
Therefore, with increasing temperature, one has to probe 
the system at increasing distances r >> £,h{T) (or 
L^^ « q « l/£^h{T)) in order to measure the asymp- 
totic form of the correlation function IjlUI) . This might 
hinder accurate determination of the universal jump. 



C. Discussion 

Before concluding this paper we would like to make 
a few additional remarks. We considered a pair of in- 
terfering quasi-condensates, however most of our argu- 
ments can be generalized to the case of several inde- 
pendent condensates. Of particular interest is a peri- 
odic array of tubes^aii or pancakes created by an opti- 
cal potentia&2ii2S. The interference pattern in this case 
shows correlations at a set of wavectors Qn = nQ, where 
n is an integer and Q is determined by the distance be- 
tween neighboring condensates. The size and angle de- 
pendence of the average interference amplitudes for each 
of these wavevectors should have the same scaling proper- 
ties as two quasi-condensates. However, the distribution 
function of fringe amplitudes will be different. In par- 
ticular, in the limit of a large number of condensates, 
the distribution function should become very narrow. 
This follows immediately from the observation that in 
this limit, higher order correlation functions in TOF im- 
ages are dominated by products of two point correlation 
function in different condensates, so there should be no 
broadening associated with Eq. H15(l below. 

Another point worth making regards the possibility of 
making analogous experiments with Fermions. For exam- 
ple, one can consider an interference of two independent 
ID fermionic systems. One obvious difference from the 
bosonic case will be the change of sign in the correlation 
function Q , reflecting different statistics of the fermions 
(this corresponds to fcrmion antibunching) . More impor- 
tantly the correlation function decays as l/\x\^^'^''-^'^^^^\ 
i.e. as 1/x or faster. This means that the integral in 
Q is dominated by short distances, where the Luttinger 
liquid description is not sufhcient, and that the integral 
converges as L — s- cx). Infrared convergence of Q im- 



plies trivial scaling \Aq\ oc VL- Moreover, the inte- 
grals appearing in all moments of the distribution are 
similarly infrared convergent, which results in a Pois- 
sonian fringe distribution at large L: P{x) oc e~^, as 
found for bosons at high temperature. In order to ex- 
tract information on the Luttinger parameter one can 
analyze the decay of density-density correlations in the 
noise {pint{x, zo)pi„t{x, zq + z)) as a function of z, which 
are directly related to the integrand of Q . We note that 
these correlations have an oscillating component, similar 
to Friedel oscillations, with wavevector 2k j. The oscil- 
lating component appears as a peak (cusp singularity) in 
the angular dependence of the interference contrast at an 
imaging angle feg tan0 = 2fc/ (see Eq. Q). The shape of 
this cusp as well as of the cusp at 6* = holds information 
on the Luttinger parameter. A more detailed analysis of 
the fermionic case lies beyond the scope of this work and 
will be addressed in future work. 

In conclusion, we analyzed the interference between 
two independent quasi-condensates. We showed that 
scaling properties of interference fringes directly probe 
the algebraic off-diagonal correlations. In particular, for 
one dimensional condensates the scaling with imaging 
length or with temperature allows to extract the Lut- 
tinger parameter. In the case of two dimensional con- 
densates this method provides a unique probe of the 
Kosterlitz-Thouless transition. We also argued that in 
the ID case one can use the distribution function of the 
interference amplitude (which is also equivalent to the 
full counting statistics of interfering bosons) as the qual- 
itative probe of the Luttinger constant. In particular, 
at if 3> 1 the distribution is narrow and at X ^ 1 or 
at finite temperatures it becomes wide Poissonian (see 
Fig. 13) . In the 2D case we expect a sharp change in the 
shape of the distribution function at the KT transition. 
The scaling analysis remains intact if more than two in- 
dependent condensates are present, but the distribution 
functions can no longer be used as a probe of the corre- 
lations. 



III. METHODS 
A. Luttinger liquid parameter 

The Luttinger liquid provides a universal long- 
wavelength description of one dimensional interacting 
Bose liquids which allows to calculate the long distance 
behavior of correlations such as (O . In certain regimes it 
is possible to derive the Luttinger parameter K and the 
healing length S^h from the microscopic interactions. In 
particular, for bosons with weak contact interactions, rel- 
evant for ultra cold atom systems discussed in this work, 
one can use Bogoliubov theory to obtainSSiSi: 
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Here 7 = 2/{aidp) < 1 is a dimensionless measure of 
the interaction strength and aid is the one dimensional 
scattering length. Analytic expressions for these param- 
eters may also be derived in the limit of strong contact 
interaction 7^1 — : 



where the thermal correlation length is 



1 



6 ^ Krnej T 



(18) 



1 + 1. 
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(14) 



B. Moments of the fringe amplitude 



Eq. H17() is valid for sufficiently low temperatures so that 
> 61, or equivalently T < h^/m(f^. For z < Eq. 
()17|l reduces to the zero temperature correlation (jSJ. In 
the opposite limit z the correlation function 117|) 

may be approximated by 



All the moments in the distribution of \ Aq\'^ can be ob- 
tained by generalizing the two point correlation function 
in Eq. JJI to the 2n-point correlation function 

(I^qI'"> - / ■•• / dzi...dz^dz[...dz'^ 
Jo Jo 

|(at(zi)...at(z„)a(z;)...a(z;))|2(15) 

For bosons with repulsive interactions described by the 
Luttinger parameter K we have 



{aHz)a{0))^p 



-Trz/2K^T 



(19) 



As we already noted, for sufficiently low temperatures 
when > L, the fringe amplitude may be found from 
Eq. 10. For i 3> equation H19|) implies 



(20) 



{\Aq\ 



2n\ 



/"^ dwi duj' Finally, substituting lfTH|) for gives ©. 

<J <J olo/~i ni-~vfd fViof ^-Vi^ o-mTnloT- Ac^•r\/r^■nA^:l•r^r 



n 



i<j 



in 



1/K 



(16) 



We also note that the angular dependence of the fringe 
amplitude at finite temperature is given by: 



which is precisely of the form (I 1 2") = (|^Q|2)"jF„(K). 
Integrals of the type appearing in Eq. IjlGI) have been dis- 
cussed by Fendley et.ali^^, who demonstrated that they 
can be calculated using special properties of Jack poly- 
nomials. From the knowledge of all moments one can, in 
principle, construct the full distribution of the interfer- 
ence fringes amplitude. In this paper we only discuss the 
limits of weak {K ^ 1) and strong [K k, 1) interactions. 



2TA/KA--k 



{\AQ\\e))^p'Le^'^i 



K/tt 



l + if2^2^2/7r2' 



(21) 



From this expression it is harder to extract K directly 
because of uncertainty in determination of and hence 



C. Finite temperature correlations in one 
dimension 
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{a\z)a{Q)) ^ pCh 



/2K ( T^lil 



\ sinh TTZ 



l/2if 



(17) 
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